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Abstract 

Speed ensemble of bistable (combustion) fronts in mean zero sta- 
tionary Gaussian shear flows inside two and three dimensional chan- 
nels is studied with a min-max variational principle. In the small 
root mean square regime of shear flows, a new class of multi-scale test 
functions are found to yield speed asymptotics. The quadratic speed 
enhancement law holds with probability arbitrarily close to one under 
the almost sure continuity (dimension two) and mean square Holder 
regularity (dimension three) of the shear flows. Remarks are made on 
the conditions for the linear growth of front speed expectation in the 
large root mean square regime. 
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1 Introduction 



We consider propagation speeds of reaction-diffusion fronts in random shear 
flows. The model equation is: 

u t = A x u + 5B -V x u + f{u), (1.1) 

where x = (xi,x) 6 D = R 1 x Q, Q = [0,L] n ~ l , n = 2,3; 5 > a scaling 
parameter; the vector field B = (b(x,u),0), b(x,ui) is a scalar stationary 
Gaussian field with zero ensemble mean and almost surely continuous sample 
paths; / is either one of the following: 

(1) the combustion nonlinearity with ignition temperature cutoff: f(u) = 0, 
11 6[0,/ 1 U6(0,l);/(n)>0 ) i t e(^l)J(l)=0 1 /GC 2 ([0,l]; 

(2) the bistable nonlinearity: f(u) = u(l — u)(u — //), fi G (0, 1/2). 

The boundary condition is zero Neumann at R 1 x dQ. The initial data Uq 
belongs to the set I s . The set I s for bistable / consists of bounded continuous 
functions with limits one and zero at X\ ~ ±oo respectively. For combustion 
/, one requires also that u decay to zero exponentially at — oo. For each 
realization of B, there exists a traveling front solution of the form u = <3>(xi + 
ct, x), with unique speed c and profile $ jlj. For u G I s , u(x, t) converges to 
a traveling front at large times [Tlj . 

We are interested in the asymptotic behavior of c = c(S) when 5 is small 
or large. Such speed asymptotics have been recently studied for both de- 
terministic fronts El El El El EU E] and KPP fronts in random shears 
[T2l IT7j . See also [12] for related applied science literature. The bistable 
(combustion) front speed c has a min-max variational characterization jTUHJ- 
Following the notation of [8 , define the functional: 

/ / \ ,/ / » _ Lv+ f(v) _ Av + 6b{x)d xi v + f{v) 
4,{v) = Mv(x)) = ^ = . (1.2) 

The strong form of the min-max variational formula of |H] for c(8) is: 

sup inf i/j(v(x)) = c(S) = inf supi[)(v(x)) (1.3) 

veK x€D v£K xeD 

where K is the set of admissible functions: 

K = {v e C 2 (D) \ d Xl v > 0, < v(x) <l,v e I.}. 
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The variational principle (jl.3|) is a powerful tool for getting tight bounds 
on c, if one can construct an admissible test function v that approximates 
well the exact traveling front solution $. This is the case when 5 is small and 
the test function is a perturbation of the known traveling front solution when 
5 = 0. In 8j, a form of such test function is proposed for deterministic shear 
flow, and the quadratic speed enhancement law ( [13] and references therein) 
is recovered if 5 is small enough. To handle an ensemble of shear flows in the 
random case, one must know which norm of b controls the smallness of 5 by 
analyzing an admissible test function v and the functional ip(v). A new form 
of multi-scale test function is found for this purpose, and it cures a delicate 
divergence problem when perturbing the traveling front at 5 = 0. It turns 
out that the smallness of 5 for the quadratic speed enhancement depends on 
the maximum (Holder) norm of b if n = 2 (3). 

Our main results are that the quadratic speed enhancement law holds 
with probability arbitrarily close to one, under the almost sure continuity 
(n = 2) or mean square Holder regularity (n = 3) conditions of the shear. As 
b can be arbitrarily large (though with small probability), our admissible test 
functions are not valid for the entire ensemble. We are unable to discuss the 
expectation of front speed. This is in contrast to the Kolmogorov-Petrovsky- 
Piskunov (KPP) case (e.g. when f(u) = u(l — u)), where there is less 
restriction on test functions and the minimal speed variational principle is 
minimization over functions on the cross section Q |12j . 

The paper is organized as follows. In section 2, the new multi-scale test 
function is presented and used to estimate the front speed and identify the 
smallness of 5 in terms of proper norms of b. In section 3, probabilistic esti- 
mates based on Borell's inequality and Karhunen-Loeve expansion of mean 
zero Gaussian fields help to control the maximum (Holder) norms of b, and 
lead to the quadratic speed enhancement laws. Remarks are made on the 
linear front speed growth law when 5 is large. 

2 Deterministic Front Speed Asymptotics 

Let us consider b to be a given Holder continuous deterministic function, 
and assume that < b >= -p- J Q b(y) dy = 0. Otherwise, the integral mean 
contributes a linear term 5 < b > to the front speed. Let U = U (x\ + Cot) be 
the traveling fronts when 5 = 0, satisfying the equation: 

U" - c U' + /([/) = 0, 
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U(—oo) = 0, U(+oo) = 1, U' > 0, U approaches zero and one at exponential 
rates. 

In the small 5 regime, define the test function 

v(x) = U{Z) + Pu((,x), (2.1) 

where the variable: 

S={l + aP)x 1 + 5 X , (2.2) 
with a a constant to be determined and x — x(%) solution of: 

-AxX = b, iGfl, 

subject to zero Neumann boundary condition at dQ. We normalize x by 
setting f n xdy = 0. The test function (|2.1|) reduces to the one in jH] if 
a = 0. However, we shall see that a ^ is essential in suppressing certain 
divergence arising in the evaluation of ip(v), similar to curing secular growth 
in multiple-scale perturbation expansions. 

For Holder continuous b(x), x is C 2+p , p G (0, 1). Assume that u is C 2 , 
and decays sufficiently fast at £ infinities. A straightforward computation 
shows that: 

Lv + f(v) = U" + f(U) 

+6 2 (2 a U" + A^u + U"\V X \ 2 + f(U)u) 

+5 3 (baU' + 2V 5 %- Vx) 

+5 4 (a 2 U" + 2 a % + A 5 u \ V X \ 2 ) 

+ f( v )-f(U)-f(U)5 2 u 

+5 5 abu^ + 5 6 a 2 u^, (2.3) 

where f(v) — f(U) — f'(U)5 2 u = 0(5 4 u 2 ) and the 0(5) terms cancel by choice 
of — Ax = b. On the other hand, up to an undetermined constant 7, we have: 

(c + 5 2 ^) d Xl v = c U' + 5 2 (c aU' + c U£ + -ylJ') 
+5 4 ('jaU' + Cocrug + 7^) 
+5 6 (a-fu ( ). (2.4) 

All terms in (|2.3|) and ()2.4|) are evaluated at (£,x). Let us choose u to solve 
the equation: 

Lu = A^u-cou^ + f (U)u = -\Vx?U" -2aU" + 1 U' + c aU', (2.5) 
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subject to zero Neumann boundary condition at dQ, and exponential decay 
at £ infinities. The solvability of (|2.5|) will be discussed later. Then we would 
have 

Lv + f(v) = (d Xl v)(c + 5 2 7 ) + R x . (2.6) 
The remainder R\ is: 

Rx = S 3 A + S 4 B + S 5 D + S 6 E + 0{8 4 u 2 ), (2.7) 

where 

A = ball' + 2V 5 m c • Vx 

B = a 2 U" + 2a + u | Vx| 2 — 7 & U' — cq a u% — 7 
D = abu^ 

E = a 2 u^ — a^u^. (2.8) 

The constant 7 is determined by a solvability condition for ()2.5j) . The 
right hand side must be orthogonal to the function U'(£) e~ c °^ which spans 
the kernel of the adjoint operator L*v = A^v + c v^ + f'(U)v. Thus, 

7 / {U'fe-^di = -L f ((|Vx| 2 + 2a)U"U' - c a(U') 2 ) e^d^dx 

Jr |"| Jd 

= ^(|V X | 2 + 2a) I (U'fe-^dt-coa I (U'fe^ d£ 
1 Jr J r, 

where the braket denotes the integral average over Q. We have performed 
integration by parts once, the boundary terms decay to zero and the integral 
f R (U') 2 e~ c °£ d£ converges for both bistable and combustion type nonlineari- 
ties. It follows that equation (J2.5|) is solvable when 

7 = ^|^|Vx| 2 ^ = |(|Vx| 2 ), (2.9) 

regardless of the choice of a. For such 7, we may choose a = — |(|Vx| 2 ) so 
7 + Co a = and the w equation becomes: 

Lu = A^u - + f(U)u = (<|Vx| 2 > - |V X | 2 ) U". (2.10) 
Note that the right hand side of ()2.1())1 has zero integral average over fl. 
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The solution u of (|2.1U|) can be expanded in eigenfunctions of the Lapla- 
cian A £ : 



where Z™ -1 denotes nonnegative integer vectors in R n ~ x with at least one 
positive component; — A £ 0j = \j<pj with zero Neumann boundary condition 
on dQ. The eigenfunctions <pj are normalized so that ||0j||L°° — 1; an d 
the eigenvalues Xj = 0(|j| 2 ) for large Clearly, 4> = 1 and Ao = 0, 
= (|Vx| 2 )- The functions Uj solve the equations 

LjUj = Uj — XjUj — cqUj + f'(U)uj = —(aj — Sjfi ao) U" . (2-H) 

The j = equation has zero right hand side, and we take u = 0. If 
a were zero, equation ()2.5|) implies that the right hand side of the j = 
equation would be — a U" + 7 U' . The general solution (zeroth-mode) for uq 
is — y£ ^'(0 + const. U'(£). The linear factor £ would render v(x) unable to 
stay in the set K due to divergence at large £ (and x\). 

Now consider the equations for j 7^ 0. As shown in [T5|, the operators Lj 
are invertible in part due to Xj > for j 7^ 0. As in Lemma 2.3 of [T5] . we 
have the regularity estimate: 

Iklb(tf) < Cxlayl, jVO, (2.12) 

where the constant Ci is independent of j and depends only on U and its 
derivatives. 

Next we show that for j 7^ 0, the ratio ^7 is bounded by C2 where C2 
is a positive constant independent of j. It suffices to prove a uniform bound 
if the cij factor is one. Consider the function w j = (3U' — Uj for some positive 
constant (3. For some r > sufficiently large, f'(U) < whenever |£| > r . 
Choose /3 as 

P>Pi= ■ Cl rfW , (2-13) 

so that by (|2.12j) . u>j > in the region |£| < r . By equation (|2.11|) for Uj 
and the fact that U' solves LjU' = —XjU' < 0, we have 

Lj Wj = -fiXjU' + U". (2.14) 
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Using the fact that U" is bounded by a constant multiple of U', and that 
infj^o Aj > 0, we may further increase /3 if necessary to ensure that LjWj < 
in the region |£| > r$. Maximum principle implies that w > for all £ e i?. 
Repeating the argument for — Uj and taking into account the a,j factor give: 



u> 



<C 3 \ aj \, (2.15) 



where the constant C3 depends only on U, and (|2.15j) holds uniformly in 
j 7^ 0. Inequality (|2.15j) says that there is no resonance when inverting the 
operator Lj to find Uj, for j 7^ 0, in other words Uj decays same as U' at 
infinities. 

To improve the estimate above, we move the term f'(U) Uj to the right 
hand side of (|2.11|) which is then bounded by C' 3 \cij\U', for C 3 independent 
of j. The left hand side operator becomes it'- — XjUj — c it'-. An explicit 
formula can be written for Uj, and the asymptotics of Xj imply that 



v® - Ci (Trm- (2 - 16) 

for a constant C4 depending only on U, where i = 0, 1,2 denotes the order 
of £ derivatives. 

Suppose that b is Holder continuous, then Schauder estimates give 

ll|Vxri| a+p <C||C (2-17) 

for some p G (0,1), C = C(Q), \\ ■ \\ p is the standard Holder norm. For 
rectangular cross section Q (dimension n — 1), the tfi/s are trigonometric 
functions. Then 

iflii<CBii6ii;(i+ijirt i+p J, p g(o,i), (2.18) 

for some constant C 5 jTHJ - Combining (|2.16j) . (|2.17j) . and (|2.18jl . we see that 
the eigenfunction expansion of u 

fi(e,5)= E "AOoM) (2-19) 
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converges uniformly in (£,£). Moreover, 



< E i^mi^m/u' 



< c 4 c 5 \\b\\i e (i + bir (3+p) = c 6 



(2.20) 



if 71 = 2, 3. The mixed derivative term is bounded as: 

|V ^ }I < E l ^\mm 



< C,\\b% E (l + bl)" (2+p) = C 7 ||fe||?, (2.21) 



if n = 2, 3. In case n = 2, f2 is an interval, then 6 e L°°(f2) suffices, because 
x is Lipschitz, so |a 3 -| < 0(||6||^ |j| _p ) for some p 6 (0, 1) which replaces 
estimate (j2.18|) . The exponent in (j2.21|) goes down to 1 + p, yet enough for 
convergence on j G Z\. As a result, for n — 2, i — 0, 1, 2, the estimates: 



<G>HC, |V ^£ ,a)l <q>llC, 



(2.22) 



hold. 

As a consequence of these estimates, let us verify the admissibility of 
test function v. Clearly, v approaches zero (one) exponentially at X\ = — 00 
(+00). As 

v Xl = (1 + 5 2 a)U' + 5 2 {l + 5 2 a)u^ 

we have for n = 3: 



2 2 " " p 



U'>-U'> 0, 



if: 



5 2 <min(^,i(C 6 ||6||5) 



2\-l 



For n = 2, in view of (J22D, > \u' if 



5 2 < min 



1 1 

2kd' 2 



(C 8 



12 \-i 



(2.23) 
(2.24) 

(2.25) 



It follows from the monotonicity of v in x\ and its asymptotics near X\ 
infinities that < v < 1 for all Xi,x. 

With v being an admissible test function, we see from ()2.6|) that 

^WHr^M^ + ^ + A-, (2 . 26) 

C/jjj (7 U Xl V 

with i2j defined by (j2~7|) . By (l2~2~l) or (l2~23J) . > Ji7'. Also, each term 
in i?i is bounded by a multiple of U'. It follows that for S satisfying ()2.24|) or 
(Q, |#^| < |4f | < C 10 <5 3 , where C 10 = C 10 (C/,||&y,pG (0,l),ifn = 3, 
Cio = Ci (U, ||6||oo) if n = 2. We have shown: 

Proposition 2.1 Lei n = 2 or 3 ; and Zet t/ie shear flow profile b(x) be 
Holder continuous with exponent p G (0, 1) over rectangular domain Q C 
i? n_1 . There is a positive constant So depending only on U and the Holder 
(maximum) norm of b for n = 3 (n = 2) such that if 5 < 5q: 

C{5)=C0 + WT\ I \Vx\ 2 dy + 0(5 3 ). (2.27) 

I I t/ O 

A specific form of So follows from \2. 24\ ) and 

At large 5, numerical evidence suggested that lim^oo c(5)/5 exists for 
bistable and combustion nonlinearities [TT] . Such a limit holds for the KPP 
nonlinearity [Sj. 

3 Random Front Speed Asymptotics 

Let us consider b = b(x,u)) as a stationary mean zero Gaussian field with 
almost surely continuous sample paths in dimension one and two for n = 2, 3 
respectively. We are interested in the restriction of b over Q. In order to apply 
results of the previous section on each realization of 6, let us write b = b + bi, 
where b = f n b(x,u) dx. Correspondingly, c = c(5,u) = co — Sb + c 2 (uj), 

where c 2 depends on b\. As the size of So depends on either the maximum 
norm or the Holder norm of b± which is an unbounded random variable, there 
is no uniform way to choose a 5q for all realizations. On the other hand, the 
probability of the occurance of very large maximum (Holder) norm of &i 
is often small, so the random speed asymptotics may hold with probibility 
arbitrarily close to one. Let us consider n = 2 first. 
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Theorem 3.1 (Two dimensional channel) Let D = R 1 x [0,L], = 
[0, L\, and b(x) be the restriction on Q of a mean zero stationary, Gaus- 
sian random process with almost surely continuous sample paths. Let f be 
a bistable or combustion nonlinearity. Then for each small e G (0,1/4) 
and q G (0,1), there is a deterministic constant 5 = <5o(e, l) such that if 
Se(0,5 ), 



Prob 



c(6,u)-c + 5b-6 2 -(( / b 1 (y)dy) > 



> K5 3 - q \ < e. 



where the constant k > is independent of 5, e, and q. As e — * 0, the 
constant 5q can be chosen to satisfy 



ao(e,g)>C|log(e)n 

Proof: Recall the inequality (Lemma 3.1, p. 62, j2]): 

E[sup b(x)} < Esup \b(x)\ < E\b(0) \ + 2E[sup b{x) 
n n n 



(3.1) 



(3.2) 



and the Borell inequality on mean zero Gaussian process with almost surely 
continuous sample paths (Theorem 2.1, p. 43, [2]): 



H = E[sup b(x)} < oo 
n 

P(|sup b(x) -fj,\ > A) < 2e 



-A 2 /(2(T 2 ) 



(3.3) 



where a 2 = E[b 2 }. With the choice of A = a/-2ct 2 log(e/4), inequality (|33|) 
implies that 



> A + fi) < P(sup b(x) - fi > A) + P(- inf b(x) — // > A) 
n n 

< 4e" A2 /( 2CT2 ) < e. 



(3.4) 



We now see that for e > sufficiently small, A > 1 and there is a constant 
Cu independent of e such that if 



5 < 5 (e) = C u 



H + A 



< C 



ii 



1 

H + A 



(3.5) 
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then (|2.25|) holds with probability at least 1 — e. Therefore, by Proposition 
2.1, the front speed asymptotics hold with probability 1 — e. Using (|2.7|) . 
()2.8|) . and ()2.22|) we see that the remainder in the expansion satisfies 

<C a S s \\b\\i, < C 12 5 3 ^ + \r 

< C 13 5 3 -^(/i + A) 4 

< C 13 5 3 ~ q 5% (n + A) 4 

< C 14 5 3 - q 

with probability 1 — e. The estimate on So(e) follows from ()3.5|) . 

Remark 3.1 An example of such b process is the mean zero stationary 
Ornstein- Uhlenbeck (O-U) process which is both Gaussian and Markov. The 
O- U sample paths are almost surely continuous. 

The n = 3 case requires a probabilistic estimate of Holder norm of b. For 
a mean zero Gaussian field, this is related to the structure of the covariance 
function R(t, s) = E[b(t)b(s)], t, s 6 Q. If the covariance function is continu- 
ous, positive and non-negative definite, there exists a Gaussian process with 
this covariance [2]. The symmetric integral operator: <fi — > J n R(t, s) 4>(t) dt 
generates a complete set of orthonormal eigenfunctions <pj on L 2 (Q) with 
nonnegative eigenvalues Xj, j = 1, 2, • • • . Define: 

p(u)= max [E\b(s)-b(t)\ 2 ] 1/2 . 

\\s-t\\<\u\V2 

Consider the partial sum (m a positive integer): 

m 

X^ m \t,u) = y/XjMWii"), ( 3 - 6 ) 
i=i 

where fy's are independent unit Gaussian random variables. The convergence 
of the partial sum to b is given by Garsia's Theorem (Theorem 3.3.2, p. 
52, 1J). It says that if (— logu) 1//2 dp{u) < oo, then with probability 
one, X^ m \t) are almost surely equicontinuous and converge uniformly on 
Q. The resulting infinite series is the celebrated Karhunen-Loeve expansion. 
Moreover, the following estimate holds for all m: 

\X^ m \s) - X^ m \t)\ < l6V2[\ogB\y 2 P (\\s-t\\) 

r\\s-t\\ 

+32^2 {-\ogu) 1/2 dp{u), (3.7) 

Jo 



Ri 



d xi v 
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where B = B(u) is a positive random variable, E[B 2 ] < 32. 

Suppose that p(u) is Holder continuous with exponent s G (0, 1), then b 
is almost surely Holder continuous with exponent s, and the Holder norm of 
b is bounded by ai[logi?] 1//2 + ct2, where oti, cti are two positive deterministic 
constants. Chebyshev's inequality gives: 

Prob ([log 5] 1/2 > A) = Prob (B > e X *] 
This implies that 



< E(B 2 )/e 2 ^ 2 < 32 e 



-2A 2 



Prob (||6||, > A) < 32e ~'\—) < C 15 e~ 2X <e 



(3.8) 
(3.9) 



if A = ^/log(^-). For e > sufficiently small, A > 1 and we take 



5 < <We) = C 



< C 



so that ()2.24j) holds with probability at least 1 — e. As before we use (12.7)1 . 
(12. 8j) . and ()2.22|1 to conclude that the remainder in the expansion satisfies 

Ri 



d Xl v 



< Crf I 



\t < C 16 5 3 - q 



with probability 1 — e. We now conclude from Proposition 2.1: 

Theorem 3.2 (Three dimensional channel) Let D = i? 1 xfi, f2 = [0, L} 2 , 
and b be a mean zero Gaussian process such that the function p(u) is Holder 
continuous. Then b has almost surely Holder continuous sample paths. For 
each small e G (0,1/5) and q G (0,1) there is a deterministic constant 
5o = 5o(e, q) such that if 5 G (0, So), 

1 



Prob 



c(6,u) -c + 5b-5 2 -{\V iX \ ) 



> k8 



3-g 



< e, 



where x satisfies: —A^x — b\, subject to zero Neumann boundary condition 
at dfl. The constant k > is independent of 5, e, and q. As e — > ; the 
constant 5q can be chosen to satisfy 



5 (e,q) >C|log(e) 



(3.10) 



Remark 3.2 Suppose that lim^oo c(5,u)/S exists almost surely and that 
^(IHloo) is finite. The dominated convergence theorem implies that 
lim^oo E[c]/5 exists and is finite. This argument is same as for the KPP 
case, see J7IJ/ for details. 
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4 Concluding Remarks 



Bistable and combustion front speeds in mean zero Gaussian random shear 
flows have been studied with the min-max variational principle of [H|. The 
quadratic speed enhancement law is valid with probability arbitrarily close 
to one in both two and three dimensional channels under the almost sure 
continuity and the mean square Holder regularity conditions of the Gaussian 
shear flows. 

It would be interesting to extend results here to the case of general convex 
cross section Q with smooth boundary, or to the case of Q in dimension higher 
than two. 
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